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We derive an ito's-type formula for the one dimensional stochastic heat equation 
driven by a space-time white noise. The proof is based on elementary properties of the 
5-transform and on the explicit representation of the solution process. We also discuss 
the relationship with other versions of this Ito's-type formula existing in literature. 
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1 Introduction 

Consider the following stochastic partial differential equation (SPDE): 



I u t (0) = u t (l) = 0, u = v ' ' 

where d t := f, d xx := Jj, W t;X : = ^gf and {W t)X ,t G [0,T],x G [0,1]} is a Brownian 
sheet. By solution to this equation we mean an adapted two parameter stochastic process 
{u t (x),t e[Q,T\,xe [0, 1]} such that for t G [0,T] and x G [0,1], 

u t {0) = u t {l) = 0, u (x) = 0, 

and such that for all / G Cq(]0, 1[) the following equality 

1 r l ft 



holds for t G [0,T]. Here (, ) denotes the inner product in £ 2 ([0, 1]) and 
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It is well known (see e.g. [8]) that equation (1.1) has a unique solution {u t (x), t G [0, T], x G 
[0, 1]} which is continuous in the variables (t, x) and that it can be represented as 



where {gt(x, y), t G [0, T], x, y G [0, 1]} is the fundamental solution of the heat equation with 
homogenous Dirichlet boundary conditions, i.e. 



Since for fixed x G [0, 1] the process t ^ u t (x) is not a semimartingale, we can not deal with 
it by means of the classical stochastic calculus. It is therefore natural to ask whether an 
Ito's-type formula can be found for this kind of process. Two recent papers, [2] and [9], are 
devoted to the investigation of this problem. In [2] the authors develop a Malliavin calculus 
for the solution process u t (x) in order to obtain an Ito's-type formula whose proof makes 
also use of projections on Wiener chaoses of different orders. In [9] the author applies the 
ordinary Ito's formula to a regularized version of the solution process u t (x); then he studies 
the limit when that regularized process converges to the real one. The main feature of this 
procedure is the appearance of a renormalization of the square of an infinite dimensional 
stochastic distribution. 

The aim of the present paper is to propose an alternative approach (and a corresponding ver- 
sion of the Ito's-type formula) to the above mentioned problem which is somehow in between 
the techniques of the articles [2] and [9]; in fact we utilize notions of white noise analysis, 
an infinite dimensional stochastic distributions theory, and prove the resulting formula via 
scalar products with test functions, more precisely stochastic exponentials. The key point is 
the gaussianity of the solution process and its explicit representation as a stochastic convo- 
lution. Our formula is in the spirit of the Ito's-type formula for gaussian processes derived 
in the paper [7j. Moreover the idea of using the properties of the semigroup associated to 
the one dimensional Brownian motion (see the proof of Theorem 2.2) is analogous to the 
one used in [Sj. 

The paper is structured as follows: in Section 2.1 we recall basic notions and facts from the 
white noise theory; then in Section 2.2 we state and prove our main result; finally in the 
concluding section we compare our formula with those already existing in the literature. 

2 Main result 
2.1 Preliminaries 

In this section we fix the notation and recall some basic results from the white noise theory. 
For additional information about this topic we refer the reader to the books [3] and pE] and 
to the paper pQ. 

Let (n,^ 7 , V) be a complete probability space and {W tjX ,t G [0, T], x G [0,1]} a Brownian 
sheet defined on it. Assume that 




d t gt(x,y) = \d xx g t (x,y) 

9t(0, y) = #t(l, y) = 0, g Q (x, y) = 5(x - y) 



(1.2) 



F = a(W t , x ,t G [0,T],x G [0,1]) 
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so that each element X e £ 2 (Q, J 7 , V) can be decomposed as a sum of multiple Ito integrals 
w.r.t. the Brownian sheet W, i.e. 

X = I n (h n ) (convergence in £ 2 (f2, JF, P)), 

n>0 

where for n > 0, /i n is a deterministic function belonging to £ 2 (([0,T] x [0, l]) n ) and I n (h n ) 
is the n-th order multiple Ito integral of h n . 

For example, the solution of the SPDE (1.1) has the decomposition: 

u t (x)= / g t - s {x,y)dW S! y = h(l m (-)g t -.(x, ■)). 
Jo Jo 

If A : D(A) C £ 2 ([0,T] x [0, 1]) -> £ 2 ([0,T] x [0, 1]) is the unbounded operator 



Ah(t,x) := A t A x h(t,x) := \[-Qp\J -Q^ h ^ x ) 
with periodic boundary conditions, we define its second quantization as the operator, 
F(A) : D{F(A)) d£ 2 (n,F,V) -> £ 2 (Q,F,V) 



^ T{A)(Y,In(K)):=Y,In{A® n h n ). 

n>0 n>0 n>0 

For p > 1 the space 

(5) p := {X e £ 2 (Q,F,V) s.t. E[\T(A p )X\ 2 } < +00} 

is a subset of £ 2 (Q, J 7 , V) and if q > p then (S) g C (5%. The Hida's test functions space is 
defined as 

(S) := f|(5) p . 

P >i 

Its dual w.r.t. the inner product of £ 2 (f2,jF, P) is called Hida's distributions space and 
denoted by (£)*. It can be shown that 

(sy = U(s)- p , 
p>i 

moreover by construction 

(s) c c 2 (n,F,v) c (sy. 

For / G C °°(]0,T[x]0, 1[), the random variable 

£ T (f) := exp{ f f f(s, y)dW s , y - \ I f f 2 (s, y)dyds}, 
Jo Jo A Jo Jo 

belongs to [S)\ therefore if ((,)) denotes the dual pairing between (S)* and (S) then the 
application 

Xe(Sy^S(X)(f):=((X,£ T (f))) 
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is well defined and it is called S-transform. The quantity S(X)(f), f G Cq°(]0, T[x]0, 1[) 
identifies uniquely the Hida's distribution X. In particular, given X, Y G (>?)*, we denote 
by XoF the unique element of (S 1 )* such that 

5(IoF)(/)=5(I)(/)5(7)(/), for all / G C o °°(]0, T[x]0, 1[); 

the stochastic distribution X o Y is named McA; product of X and F. 
Notice also that if X G C 2 {VL, T ,V) then 

S(X)(f) = ((X,£ T (f))) = E[X£ T (f)}, 

and if £ is an It 6 integrable stochastic process then 

5( / T f\t, x dW ttX )(f) = E\ f [\ t , x dW t , x £ T (f)} = f f E[^ x £ T (f)]dxdt. 
Jo Jo Jo Jo Jo Jo 

2.2 Ito's-type formula 

Before the main theorem of this paper, we state and prove the following auxiliary result. 
Lemma 2.1 Let u t (x) be the solution of the SPDE (1.1); then d xx u t (x) G (S)* . 

Proof. It is known that the fundamental solution g t (x,y) of the heat equation with 
homogenous Dirichlet boundary conditions (1.2) can be represented as 

9t(x,y) = 5^V A "*e n (x)e n (|/), 

n>l 

where e n (x) = \f2sin{imx) and \ n = ir 2 n 2 . This gives 

d xx g t (x,y) = ^2e~ Xnt (-\ n )e n (x)e n (y) 

n>l 

= ^2e- Xnt e n (x)d yy e n (y) 

n>l 

= J2 e ~ Xnte n(x)(-A 2 y e n (y)). 



n>l 



Therefore, 

E[\T(A- 2 )d xx u t (x)\ 2 } = f f 1 \A; 2 A y 2 d xx g t _ s (x,y)\ 2 dyds 

Jo Jo 



t r\ 

2„ /„ „.\|2. 



JO 
t pi 



\A S g t - s (x,y)\ dyds 

I I \9t-s(x,y)\ 2 dyds < +oo, 
Jo Jo 

proving that d xx u t (x) G (S)- 2 C (S)*. □ 



< 
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Theorem 2.2 For any if G C 2 (R) and I G CqQO, 1[) one has 

(tp(ut),l) = (<p(u ),l) + J (ip'(u s )l,dW s ) + \(J o v'(u s )od xx u s ds,l) 



+ l([\"(u s )da 2 (s;.)J), (2.1) 



2 

where 

• a 2 (s,x) = E[\u s {x)\ 2 ] = f* J* g 2 s _ v (x,y)dydv; 

• J Q * Lp"(u s (x))da 2 (s; x) is a Stiltjes integral in the variable s; 

• ip'(u s ) od xx u s is the Wick product between (p'(u s ) and d xx u s . 

Proof. We aim at proving that 

EMu t ),i)s T (f)] = E[ns T (f)], 

for all / G Co°(]0, T[x]0, 1[) where 1Z denotes the r.h.s. of equation (2.1). This fact together 
with the property 

span{£ T (f),f G C o °°(]0,T[x]0,l[)} is dense in C 2 (n,F,V) 

will imply that 

{ip(u t )i I) — 1Z V — almost surely, 
which is the statement of the theorem. 

Let us fix an arbitrary / G C£°(]0, T[x]0, 1[); a simple application of the Girsanov's theorem 
yields: 

t rl 



E[ip{u t {x))£ T {f)\ = E [v(J J g t - s {x,y)dW s , y y T {f) 



E 



t r l ft rl 



/ / 9t-s(x,y)dW S! y+ / / g t - s (x,y)f(s,y)dyds 
o Jo Jo Jo 



Now observe that 

i-t rl rt rl 



/ / 9t-s(x,y)dW 3 ,y+ / / 9t-s(x,y)f(s,y)dyds, 
Jo Jo Jo Jo 

is a Gaussian random variable with mean given by 

m(t,x):= / gt- s (x,y)f(s,y)dyds, 
Jo Jo 



and variance equal to 



a 2 (t,x):=[ [ g 2 _ s (x,y)dyds. 
Jo Jo 
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Therefore we can write 



E[<p(u t (x))£ T (f)} = E ip(^ / / g t -s(x,y)dW s , y + / g t - s (x,y)f(s,y)dyds^J 



'0 JO JO Jo 

(P a 2 {t , x) <p)(m(t,x)), (2.2) 



where {Pt}t>o denotes the one dimensional heat semigroup, 



(P t <p)(x) := / <p(y)-± 

JK V * 



e 2t ay. 

2ixt 



Identity (2.2) turn out to be very useful; in fact we can apply the ordinary chain rule and 

get 

E[<p{u t {x))E T {f)\ = (P a2(t , x) <p)(m(t,x)) 

f* d 

= <p(u (x))+j —(P a 2 {VjX) <f)(m(v,x))dv 
= (p(u (x)) + J (P a 2 iV)X) p') (m(v, x )) dm ^ ^ dv 
+ / ^( p a 2 {v,x)<f")( , m(v,x))da 2 (v;x) 



Recalling that 



we have 



lo 2 

=: <p(u (x)) +l(t,x) +ll(t,x). 



9v- s (x,y)f(s,y)dyds, 



f(v,x)+ [ [ \d xx g v - s (x,y)f(s,y)dyds 



dv Jo J 2 

™ pi 



f( v > x ) + \ d **{ [ J o 9v-s(x,y)f(s,y)dyds^, 



and hence 



I(t,x) = J (P a 2 iVjX) (p')(m(v,x))^f(v,x) + ^d xx (^J J g v _ 8 {x,y)f{s,y)dydsj}dv 

= / (P^(v,x)<f')(m(v,x))f(v,x)dv 
Jo 

+ J^{PaHv,x)^'){'m(v,x))]<) xx {^J^ g v _ s (x,y)f(s,y)dyds^dv 



t pt 



E[<p'(u v (x))£ T (f)]f(v,x)dv + J E[ip\u v {x))£ T {f)\^d xx E[u v {x)£ T {f)\dv, 

where the last equality is due to identities (2.2). Moreover since from Lemma 2.1 d xx u v (x) G 
(S)* we have 

d xx E[u v (x)£ T (f)] = d xx S{u v {x)){f) = S(d xx u v (x))(f), 
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where S(d xx u v (x))(f) denotes the iS-transform of d xx u v (x). 

If now I G Cq (]0, 1[) is a test function in the space variable x, by the properties of the Wick 
product we conclude that 



1 rt 



I(t, x)l(x)dx 




'0 JO 

«1 rt 



+ 




'0 JO 

-I rt 




E[(p'(u v (x))£r(f)]l(x)f{v, x)dvdx 
E[(p'(u v )£ T (f)]^S(d xx u v (x))(f)l(x)dvdx 
E[(p' (u v (x))£T{f)]l( x ) f i v '> %)dvdx 



'0 JO 

"1 /-t 



+ 




./0 



S(tp'(u v (x)) o -d xx u v (x))(f)l(x)dvdx 



E 



t rl 




tp'(u v (x))l(x)dW v - 



+ 



JO 
1 , rt 



'0 x JO 

Looking again at identities (2.2) we also discover that 

"*1 



^'(it„(x)) O -d xx u v (x)dv)l(x)dx)S T (f) 



ll(t,x) = J -E[ V "(u s (x))S T (f)]da 2 (s;xy, 

combining the expressions for X and XX we can now conclude that 

(E[ip(u t )&r(f)],l) = (<p(u ),l) 
+ E 




+ E 
+ E 

This completes the proof. 



' rl f* 1 

V5 / (w 1 ,(x))/(x)<iW / 1 , i:r + ( / <p'(u v ) o -d xx u v dv, I) )S T (f) 

o Jo Jo * 



+ (J -E[(p"(u s )£ T (f))da 2 (v; •), I) 

= (ip(u ),l) 

t ,t x 

(tp'(u v )l, dW v ) + ( / o -d xx u v dv, I) )S T {f) 

o Jo 

(i !\"{u s )da 2 {v--),l)]e T {f) 



□ 



2.3 Comparisons 

Zambotti's formula : In [U] the author obtains the following Ito's formula: 

((p(u t ),l) = (<p(uo),l) + ^ [ (l",tp(u s ))ds+ [ ((p'(u s )l,dW s ) 

* Jo Jo 

(I, : \d x u s \ 2 : Lp"(u s ))ds, 



(2.3) 
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where the last term is defined as the limit of renormalized diverging quantities. The 
procedure to derive this formula is to approximate the solution of the SPDE via the 
smoother process 

u\(x) •■= / / g t - s +e(x,y)dW St y, 
Jo Jo 

and then to pass to the limit for e — > 0. We are now going to show how to manipulate 
formula (2.1) to make it look like (2.3). 

Following the same line of reasoning explained in the proof of Theorem 2.2 with u\{x) 
instead of Ut(x), one can easily see that everything can be carried in a similar manner; 
the only difference consists in the fact that in this case the function 

a 2 t {t,x) := E[\u e t (x)\ 2 } = [ [ g 2 _ s+e (x,y)dyds, 

Jo Jo 

is differentiable w.r.t. t; therefore the last term of formula (2.1) can be rewritten as 
1 -J\»(ul(x))da 2 (s;x) = \ j\"(ut(x))^^-ds 
^'«(^))(^ gl{x,y)dy + J j d B (g*_ v+e (x,y))dydv)ds 



1 

2Jo 

\ rt , /■! rs r l 



v'\ u \{ x ))[ / ge( x >y) d V + g s ^ v+e (x,y)d xx g s _ v+e (x,y)dydv)ds 
o v Jo Jo Jo J 

\ I f"( u l( x )) / g 2 e{x,y)dyds + - I (p"(u e s (x)) [ [ g s - v+e (x,y)d xx g s - v+e (x,y)dydvds 
1 Jo Jo 1 Jo Jo Jo 

(f"(u e s {x)) I g 2 (x,y)dyds + - I I I D Viy Lp'(u e s {x))D Vty d xx u e s (x)dydvds 
o Jo * Jo Jo Jo 

If 1 f 1 1 f* 

(p"(u € s (x)) g 2 (x,y)dyds + - / ip' (u c s (x))d xx u e a (x) - (p'(u e s (x)) o d xx u £ s (x)ds. 
Jo z Jo 

Here we used a property of the Wick product which shows its interplay with the Hida- 
Malliavin derivative, namely 

Xo / h(s,y)dW 8 , y = X / h(s,y)dW s , y - / D s , y Xh(s,y)dyds. 
Jo Jo Jo Jo Jo Jo 

See [1] and [0] for details. Therefore equation (2.1) becomes 

(<p(ut),l) = I) + Jj<p\ul)l,dW s ) + ~(J^'(ul)<>d xx utds,l) 

+\(j\"{<) j g 2 A;y)dyds,l) 



1 



+| ( / ¥'(u e s )d xx u e s - cp'(u £ s ) o d xx u e s ds, I) 



t 



M«s),i>+ / (tp'(u%dw s ) 

Jo 

+l(f V"{<) I g 2 (;y)dyds,l} + l(l ip'(ul)d xx ulds,l). 



Moreover 

<p'(u £ s (x))d xx u e s (x) = d xx p(u e s (x)) - tp"(u £ s (x))(d x u e s (x)) 2 ; 
a substitution in the previous equation gives 

(¥>(«!),/) = (<?(«§),/) + f t (<p'(ul)l,dW s ) + ±([ t <p"(ul)([ g 2 e (;y)dy)ds,l) 

Jo A Jo Jo 

= (<p(«§),0 + f{ip\ul%dW s ) - \{ f \"{ul){{d x ulf - [ g 2 e (.,y)dy)ds,l) 
Jo A Jo Jo 

1 * 



+ 2</ o V{u%l"). 

Since this identity coincides with the expression derived in [0] before taking the limit 
for e —>■ 0, the equivalence of the two formulas is proved. 

Gradinaru-Nourdin-Tindel's formula : In [2] the authors present the following Ito's 
formula: 

<p(ut) = cpfa) + [ (ip'(u s ),8u s ) + \ f Tr(e 2sA V "(u s ))ds. (2.4) 
Jo * Jo 

The term J^(ip'(u s ), Su s ) denotes a Skorohod's type integral w.r.t. the solution process 
u t (x). If we write formally 

1 t i t f 1 i du s (x) 

5u s {x) = -d xx u s (x)ds + SW S>X and / ((p'(u s ),5u s ) = / / v?'(it s ) o — ^ — cteds, 



we get 

ft r t rl 

In 1 11.-) C> — 

ds 



((p'(u s ),8u s ) =11 <p'(u s ) o dUs } X ^ dxds 



o Jo 



t r l 



l f'( u s) o (Xd xx u s (x) + efeds 

/• I -y i-t rl 

<p'(u s ) o -d xx u s (x)dxds + (f'(u s )dW sx . 



o Jo 



JO 2 JQ J Q 



This procedure, far from being rigorous, suggests some common feature between (2.1) 
and (2.4). However the identification of the "Ito's terms" of the two formulas, namely 



jT ^"{u s {-))da 2 {s- ■) and f Tr(e 2sA ip"(u s ))ds, 
doesn't seem to be straightforward. 
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